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Abstract 

' We study the motion of a spinning test particle in Schwarzschild spacetime, 

0^ \ analyzing the Poincare map and the Lyapunov exponent. We find chaotic 

behavior for a particle with spin higher than some critical value (e.g. Scr ~ 
0.64/iM for the total angular momentum J = 4//M), where ^ and M are 
the masses of a particle and of a black hole, respectively. The inverse of the 
^ ■ Lyapunov exponent in the most chaotic case is about three orbital periods, 

(<~^ \ which suggests that chaos of a spinning particle may become important in 

. some relativistic astrophysical phenomena. 

■ The "effective potential" analysis enables us to classify the particle orbits 

\^ \ into four types as follows. When the total angular momentum J is large, 

G> . some orbits are bounded and the "effective potential" s are classified into two 

types: (Bl) one saddle point (unstable circular orbit) and one minimal point 
qh, (stable circular orbit) on the equatorial plane exist for small spin; and (B2) 

two saddle points bifurcate from the equatorial plane and one minimal point 
^j^l remains on the equatorial plane for large spin. When J is small, no bound 

orbits exist and the potentials are classified into another two types: (Ul) no 
extremal point is found for small spin; and (U2) one saddle point appears on 
the equatorial plane, which is unstable in the direction perpendicular to the 
equatorial plane, for large spin. The types (Bl) and (Ul) are the same as 
I those for a spinless particle, but the potentials (B2) and (U2) are new types 

caused by spin-orbit coupling. The chaotic behavior is found only in the type 
(B2) potential. The "heteroclinic orbit", which could cause chaos, is also 
observed in type (B2). 
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I. INTRODUCTION 



Chaos is now one of the most important ideas used to explain various non-hnear phe- 
nomena in nature. Since the research on the three body problem by Poincare, many studies 
about chaos in celestial mechanics and astrophysics have been done and revealecl the im- 
portant role of chaos in the Universe |T|, 0. Although we know many features of chaos 
in Newtonian dynamics, we do not know, so far, so much about those in general relativity. 
If gravity is strong, e.g., a close binary system or a particle near a black hole, we have to 
use Einstein's theory of gravitation. Because the gravitational field in general relativity is 
non-linear, we may find a new type of chaotic behavior in strong gravitational fields, which 
do not appear in Newtonian dynamics 0- 0. In a previous paper 0, we studied a cri- 
terion for chaos of a test particle motion around an A^-black hole system (or an A^-naked 
singularity system) and found that a local instability determined by the Riemann curvature 
tensor provides us a sufficient test for chaos. We also found that the existence of an unstable 
circular orbit, which guarantees the existence of a homoclinic or heteroclinic orbit, plays a 
crucial role for chaos. However, the relativistic systems analyzed so far by several authors 



13|, in which chaotic behavior of a test particle is found, are rather unrealistic |B[- [11 



except for the perturbed spacetimes of the Schwarzschild black hole solution [|T^, |[T^. As 
for the other interesting cases of systems, for example, the A^-extreme black hole system |^- 
is unstable, the existence of a strong uniform magnetic field around a black hole |T0] is 
not likely, and naked singularities may not exist. We may wonder whether any realistic 
relativistic system can be chaotic and when chaos may play an important role in such a 
relativistic astrophysical phenomena. 

In astrophysics, rotation of a system plays a quite important role. The angular momen- 
tum or spin may completely change the evolution of the system. In a dynamical system, 
rotation or spin is one of the most important elements and it may sometime cause chaotic 
behavior. For example, the Henon-Heiles system, which describes the motion of a star mov- 



ing in the potential of a rotating galaxy, is chaotic [|1J]. We also know that some spin-orbit 
interaction induces chaos in Newtonian gravity. This may also be true in a relativistic sys- 
tem such as the evolution of a binary system. The motion of coalescing binary systems 
of neutron stars and/or black holes is very important to study because they are promising 
sources of gravitational waves, which we are planning to detect by large-scale laser interfer- 
ometric gravitational observatories, such as US LIGO HTSf. If we will detect the signal of 



gravitational waves emitted from these systems and compare it with theoretical templates, 
we may be able to determine a variety of astrophysical parameters of the sources such as 
their direction, distance, masses, spin, and so on |]T6[. In order to extract exact information 
about such sources from the observed signal we need the exact theoretical templates of the 
gravitational waveforms. To make such templates, it is very important to know the exact 
motion of sources. Hence, the equations of motion in the post- Newtonian expansion in terms 

of a small parameter e ~ {v/cY ~ GM/r have been studied by many authors Ijl^. Those 
can be written schematically as 

^"^^ - n -L I «(3/2) , „(2) , „(2) , _(5/2) ^/ (3)^ ^ 

- + apj^ + a^Q + a^p^ + ass + ^^rr J, U-J-j 

where the subscripts A^, PN, SO, 2PN, SS and RR denote Newtonian, post-Newtonian, 
spin-orbit coupling, 2nd post Newtonian, spin-spin coupling, and radiation reaction terms, 
respectively |[^; and the superscript corresponds to the order of expansion in e. To make 
sufficient templates, we may need at least the 3rd-order post Newtonian contribution to 



obtain the S/N ratio required from the observation |T^- But this is still under investi- 
gation in the world. 

The spin effect is also important. The spin terms in Eq.( p..l|) , aso and ass, induce 
a precession of the orbital plane through the spin-orbit or spin-spin coupling, resulting in 
modulation of the gravitational waveforms [^, [^]. In 0^, it is also shown that the or- 
bital plane may behave very strangely due to the spin effects. We cannot verify whether 
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or not any chaotic behavior occurs in their system. But from the studies on spin effects 
in Newtonian dynamics, we know that a spin effect can make a motion chaotic. We then 
expect that a relativistic system such as a coalescing binary pulsar may also show the similar 
non-linear phenomena. The gravitational waveform from the system with chaotic motion 
will be different from a system with a regular motion, for example, a regular precession of 
the orbital plane as shown in |^ . The chaos might be too strong to make a complete tem- 



plate of gravitational waves, or rather it might give us new information about astrophysical 
parameters from a time series of the observed waveforms. We will discuss this problem for 
a coalescing binary system with highly spinning bodies elsewhere. 

Thus we believe that a study about spin effects on the orbital evolution of a relativistic 
system and its gravitational waveform is very important from the viewpoint of observations 
as well as of academic interest. In this paper, to clarify the spin effect on the orbital motion, 
especially the spin-orbit interaction, we study the motion of a spinning test particle around 
a Schwarzschild black hole. So far, studies about a spinning test particle in relativistic 
spacetime have been done by many authors since the basic equations were derived by Papa- 
petrou and reformulated by Dixon |2^. Corinaldesi and Papapetrou already discussed 
a spinning test particle in Schwarzschild spacetime But, apart from the supplementary 
condition, from which they adopted a different equation from the present standard one, 
they presented the basic equations and discussed some terms with physical interpretations. 
They did not analyze the orbits in detail from a viewpoint of the dynamical system. Kerr 
or Kerr-Newman spacetime was also analyzed by several authors [^]. In [^, and 
PP| , the effective potential of the spinning particle is given and the spin effects on the bind- 



ing energy are discussed. In |^ and [Q, the gravitational waves produced by a spinning 
particle falling into a Kerr black hole or moving circularly around it is discussed and the 
energy emission rate from those systems is calculated. But in those papers they discussed 
only the case of the orbit in the equatorial plane or on the symmetric axis of the black hole. 
Since we are interested in chaotic motion induced by spin-orbit coupling here, we have to 
discuss the most generic situation, i.e., the orbital motion off the equatorial plane. 

This paper is organized as follows. In section 2 we shall briefly review the basic equa- 
tions, i.e., the equations of motion for a spinning test particle in relativistic spacetime, a 
supplementary condition and some constants of motion. We specify the background space- 
time to be a Schwarzschild black hole, then we write down those equations and introduce 
a sort of "effective potential" , which enables us to classify the particle behavior. In section 
3, performing numerical integrations, we show that chaos occurs for a highly spinning test 
particle. Summary and some remarks follow in section 4. 

Throughout this paper we use units c = G = 1. We define the signature of the metric as 
(-,+,+,+). 



II. BASIC EQUATIONS FOR A SPINNING TEST PARTICLE 



A. Pole-dipole Approximation 

The equations of motion of a spinning test particle in a relativistic spacetime were first 

derived by Papapetrou and then reformulated by Dixon [^. Those are a set of equa- 
tions: 

— = v^, (2.1) 

^ = -^R%XS'^'^, (2.2) 

=p'^y- -p-v^^^ (2.3) 

Dt 
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where T,v'^,p'^ and S^"^ are an affine parameter of the orbit, the 4- velocity of the particle, 
the momentum, and the spin tensor, respectively, r is chosen as the proper time of the 
particle in this paper, then v^v^ = —1. The multipole moments of the particle higher than 
mass monopole and spin dipole are ignored. It is called the pole-dipole approximation. 
We need a supplementary condition which gives a relation bet ween f ^ and p^, because 
is no longer parallel to in the present case. The r.h.s. of Eg. ( p.2|) denotes a spin-orbit 



coupling through a strong gravitational field. We adopt the following condition p5 

P^.S'^" = 0. 



(2.4) 



This condition is related to how to choose the center of mass in an extended body, and this 
choice gives a consistent condition [^]. Using (|2.4|) we can write down the relation between 
f ^ and explicitly, that is. 



2/i2A 



where 



A 



and 



N 



4A2/i 



1 + ■^j^Rap-fSS"'^ S'^^ , 



(2.5) 



(2.6) 



-1/2 



(2.7) 



is a normalization constant fixed by Vi^v" = —1. = p'^/fi is a unit vector parallel to the 
momentum p'^, where the mass of the particle n is defined by 



2 1/ 
/i = -PuP ■ 



(2.8) 



This system has several conserved quantities. Regardless of the symmetry of the background 
spacetime, it is easy to show that fi and the magnitude of spin S, defined by 



02 Q CfJ-l^ 



(2.9) 



are constants of motion |3^. If a geometry possesses some symmetry described by a Killing 
vector associated with the symmetry, we can show that 



(2.10) 



is also conserved p5 |. 



B. Spinning Particle in Schwarzschild Spacetime 

As for the background spacetime, we assume a Schwarzschild spacetime, i.e., 
ds^ = -f{r)de + f{r)-^dr^ + r^{de^ + sin^ 



(2.11) 
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where 



/(r) = l-— (2.12) 



with M being the mass of the black hole. Because the spaceti me is static and spherically 
symmetric, there are two Killing vector fields, and From ( p.lC ), we find the constants 
of motion related with those Killing vectors as 

E = -C(^t^ = -pt-^S'\ (2.13) 
J, = C(^) =p^- r(5^" - rS'^^ cot 9) sin^ 6. (2.14) 

E and Jz are interpreted as the energy of the particle and the z component of the total 
angular momentum, respectively. Because the spacetime is spherically symmetric, the x 
and y components of the total angular momentum are also conserved. Then we have two 
additional constants of motion as 

Jx = —pe sin (p — P(j, cot 6 cos 

WS'^^ sin^ ^ sin + rS^'' sin 9 cos ^ sin + rS"'^ cos 0, (2. 15) 

Jy = Pe COS 4> — P4, cot 9 sin 

+r2^^'^sin2^cos0 + r^'^" sin ^ cos ^ cos — rS*^' sin0. (2.16) 

Because the background is spherically symmetric, without loss of generality we can choose 
the z axis in the direction of total angular momentum as 

{J„Jy,Jz) = (0,0, J), (2.17) 

where J > 0. Three constraint equations ( |2.14|) ^( ^.16| ) with Eq. ( p.l?| ) are reduced to 



"^cot^, (2.18) 



S"-" = (2.19) 
r 

S^"- = -(-J+^^). (2.20) 



r V sin 9 

S^' {i = r, 9, 0) are fixed from Eq. (^ with Eqs. (^)~(g;5Up as 



S'^ = -— N + - JP<rj , (2-21) 

rpt \ sm 9 J 

-te _ 1 / ^ , JP't> 



S"' = — U^p^ + ^cot^ , (2.22) 
rpt \ r J 

S^^ = -^ijp,.-^^ + ^-^coi9]. (2.23) 



rpt \ sin 9 

The energy conservation Eq.( |2.13| ) is now 



E = -pt + ^Apl + ^n- ■ (2-24) 
PtV^ \ sm 9 J 



The procedure to give initial conditions for calculating the orbital evolution of the particle 
is as follows. First, we give constants of motion S, J and E with which the motion of the 
particle is bounded to a compact region of the spacetime. We will discuss how to do it in 
detail in the next subsection. We set the particle on the equatorial plane 6 = n/2 with r = tq 
and = 0. Next, we give the spatial components of the spin tensor S = (S*^"^, 5''^'", 5*^'^). S^'f' 
vanishes at ^ = 7r/2 from ( |2.18[ ), which means that the initial spin is perpendicular to the 



radial direction. The last parameter we need to choose initially is. 



a 



tan 



-1 



re 



(2.25) 



which determines the direction of the spin. Note that a = and vr denote the spin anti- 
parallel and parallel to the positive 2;-direction, respectively. In this case the orbit of the 
particle is always restricted to the equatorial plane and chaos never occurs. We assume 
TT /2 < a < 3tt/2, which corresponds to the case that the z component of spin points to the 
same direction as that of the total angular momentum. Otherwise, the particle cannot get 
into any relativistic region and chaos never occurs, as we will see later. 
From Eq. (O) with Eqs. (|2J8|)~ (|2]2|), we find that 



PtS 



.P sin^ a + jj?rQ 



cos a, 



•re 



PtS 



.P sin^ a + iJ?rQ 



; sm a. 



(2.26) 



Inserting Eqs. (|2.19|) and ( |2.2CI|) with these equations into Eq. (|2.24|) (at 9 = 7r/2), we get 
the quadratic equation for Uf. Solving it and inserting the result into Eq. ( ^.2(j|), we find 
the initial values for pt, S. pe,P(p and S*** are determined from Eqs. ( p.l9| ) ^ (|2.23|) and 
is given by use of Eq. (|2T 



C. Contour of Zero Meridian Momentum as an "Effective Potential" 

Since we are interested in chaotic behavior of a particle orbit, we analyze a test particle 
which does not escape to infinity and does not fall into a black hole. Therefore we have to 
choose appropriate parameters of the particle, E, J and S. If we have an effective potential, 
such a choice is easy. For example, for a spinless particle traveling around a Schwarzschild 
black hole with the orbital angular momentum L, whose motion is restricted to a plane (e.g. 
the equatorial plane 6 = 7r/2), the effective potential of the particle is given as 



V\r;L)= 1 + (2.27) 




Since the region where the particle with the energy E can move is given by V'^{r) < E"^, it 
is easy to choose the energy E such that the particle will move in a compact region. 

In the case of a spinning particle, because the spin-orbit coupling ( p.2| ) is not a potential 
force and an additional dynamical variable, i.e., the direction of the spin, exists, we cannot 
find any effective potential in the 2-dimensional r-6 plane. To know the region where the 
particle can move, however, we do not need to find an effective potential itself. Rather we 
need to know only the boundary of such a region, i.e., a curve in the r-6 plane where both 
p"^ and p^ vanish. From ( p.l9| ), S""^ also vanishes there, which means that the spin lies in the 
meridian plane and there remains no freedom of the spin direction (a = or tt). From ( p.8| ) 
with p^' = p^ = 0, we can set 

Pt = -fif^^'^ coshX, 

p^ = fir sin 9 sinh X, (2.28) 
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where X{r,6) is an unknown function. Inserting this into Eq. ( p.9| ) with Eqs. ( p.l^ ) 
( p.20| ), we find the equation for X as 



(/iV^ - S^f) sinh^ X - 2iiJr sin 9 sinh X + (J^ - S^)f + —J^ sin^ ^ = 0. 

r 

Finally, from Eq.( p.l3| ), we obtain the equations for such a curve as 

E = V^^)ir,9;J,S), 

where 



V{±)ir,e; J,S) = fx 



/2 coshX(±) + — 



MsinhX(±) / Jsin^ 



/^/^r coshX(±) y fir 



sinh X. 



(±) 



(2.29) 



(2.30) 



(2.31) 



and 



sinhX(±) = 



fiJr sin 6 

^2^2 _ ^2J 



± 



fx 



Wsin^^ UJ'-S')f+'-f J' sin' 9 



^2^2 _ g2j 



(2.32) 



The subscript (±) corresponds to the direction of the spin, i.e., if the z component of the 
spin is the same direction as the total angular momentum (S"^*" < 0, i.e., 7r/2 < a < 37r/2), 
we take V(_), while if it is the opposite {S'^^ > 0, i.e., — 7r/2 < a < it/2), V(+) should be 

applied. Note that S"^^ never changes its signature during the evolution. Imposing S = 
and 6 = tt/2, V(±) is reduced to the conventional effective potential for a spinless particle 
in a Schwarzschild black hole, i.e., Eq. (|2.27|) . The derivation also shows that the particle 
with energy E can move in the region of the r-6 plane such that V^±)(r, 6; J, S) < E'^. Then 

we shall call V{r, 6; J, S) the "effective potential" of a spinning particle in Schwarzschild 
spacetime. 

The typical shape of the potential on the equatorial plane is shown in Fig.l. From this 
figure, we can see that the particle in the potential V(_) can move closer to the event horizon 

(r = 2M) than that in the potential V(+). This is because the spin-orbit interaction is a 
repulsive force if the spin is parallel to the orbital angular momentum and it will balance 
with the gravity, so that a stable periodic orbit closer to the event horizon becomes possible 
without falling into the black hole. If the direction of the spin is opposite, we find the reverse. 
Such a spin -orbit interaction is induced through the gravitational interaction as the r.h.s 
of Eq. ( p.2|) , which breaks the integrability of the particle motion. Since we are interested 
in chaotic behavior, a large spin-orbit interaction may be more interesting. Hence, in this 
paper, we study only particle motion in a strong gravitational field, i.e., in the potential 
V(_). In what follows, we drop the subscript (— ). 

Once J and S are given, we can depict a contour map of V{r, 6; J, S), with which we find 
an allowed region where a particle with E can move. We find that the "effective potential" 
V is classified into four types depending on J and S. These are given in Fig.2. In Fig. 2a, 
the potential has one saddle and one minimal point on the equatorial plane, which appears 
when S is small compared with J. Since the spin effect is small, the shape of the potential 
is similar to that for a spinless particle. We call this potential type (Bl). The orbit in 
this case never becomes chaotic. Fig. 2b shows the potential which has two saddle points off 
the equatorial plane and one minimal point on the equatorial plane. The saddle points are 
located symmetrically on opposite sides of the equatorial plane. We call this potential type 
(B2). For fixed J, as 5* gets larger, type (Bl) changes into type (B2). This happens because 
a repulsive spin-orbit interaction is angle dependent. The "effective potential" s shown in 
both Figs. 2c and 2d have no bound region, although the latter case has a saddle point on the 
equatorial plane. The particle will eventually fall into a black hole. We call those potentials 
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type (Ul) and (U2), respectively. These potentials appear when J is small enough, i.e., the 
centrifugal force is too small to balance with the gravity. The type (Ul) potential is similar 
to that for a spinless particle. In the case of type (U2), however, the situation is slightly 
different. Because the saddle point marked by a cross in Fig. 2d is minimal in the r direction 
and maximal in the direction, a particle will gradually depart from the equatorial plane 
and fall into a black hole. There is a potential barrier on the equatorial plane from the 
repulsive spin-orbit interaction. 

We show in Fig. 3 what values of J and S belong to which types. Because types (B2) and 
(U2) never appear in the case of a spinless particle, we conclude that such strange behavior 
in those potentials is induced by a spin effect. As we will show later, we find that it is only 
for the type (B2) that the particle motion can be chaotic. Note that the value of S at the 
bottom end of the (U2) region is slightly smaller than 1/iM, i.e., at S* ~ 0.987/iM. Hence 

the type (U2) still appears with physically meaningful values of S and J, as we will discuss 
soon. 



III. NUMERICAL RESULTS 



A. Chaos in Schwarzschild black hole 



Using the Bulirsch-Stoer method |37|, we integrate the equations of motion numerically 
for various values of parameters, J and S and various initial configurations of tq and a. 
As for the dynamical variables, we solve all components of the position vector x^, of the 
momentum p^, and of the spin tensor S^^ , and use the constraint equations ( p^.4| ),( P^ ),( ^.9[ ), 
( p.l3| ) ~ (|2.16| ) to check the accuracy of our numerical integration. After 10'^ orbital periods. 



which is the end of our numerical calculation, the relative errors are smaller than 10~^^ for 
each constraint. The methods with which we analyze the chaotic behavior are the Poincare 
map and the Lyapunov exponent. To make the Poincare map, we adopt the equatorial 
plane {6 = 7r/2) as a Poincare section and plot the point (r, f^) when the particle crosses 

the Poincare section with < 0. If the motion is not chaotic, the plotted points form a 
closed curve in the 2-dimensional r-v"^ plane, because a regular orbit will move on a torus 
in the phase space and the curve is a cross section of the torus. If the orbit is chaotic, some 
of those tori will be broken and the Poincare map does not consist of a set of closed curves 
but the points will be distributed randomly in the allowed region. From the distribution 
of the points in Poincare map, we can judge whether or not the motion is chaotic. The 
Lyapunov exponent is another method to judge the occurrence of chaos and it also gives 
a naive estimation of the strength of chaos quantitatively. This denotes how fast the close 
orbits in the phase space will diverge in future. The Lyapunov exponent A is defined by 



A = lim - In 

t^oo t 



d{t) 



dio) 



(3.1) 



where d{t) is the distance at time t between two neighboring points in the phase space. If the 
orbit is chaotic, then A will converge to some positive value, which means that the distance 
will diverge exponentially with a typical time scale t\ = X~^. For a particle moving along a 
geodesic, we can calculate the Lyapunov exponent by integration of the equation of geodesic 
deviation; 



Dr2 



-R%XnPv'^, (3.2) 



where n'^ is a deviation vector. In the present case, however, because a spinning test particle 
does not move along a geodesic as discussed in section 2, we need another way to estimate 
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the Lyapunov exponent. Here, we have adopted the method developed by Sano and Sawada 
|38| to estimate the Lyapunov exponent. We prepare a time series of p(t), z{t), v^it) and 
where p and z are defined by p = rsin6' and z = rcosd, and calculate the Lyapunov 
exponent from such a set of data. The way to calculate the Lyapunov exponent is as 
follows. Let {xi\i = 0, 1,2, ■ ■ ■} denote a set of points on some orbit in the 4-dimensional 
phase space, i.e., Xi = {p{iAt), z{iAt),vP{iAt),v^{iAt)) where At is a small time interval. 
A set of "deviation vectors" {y^} at a point Xi is defined by 

{Ui} = {xk, -Xi\ < \\xk^ -Xi\\ < e}, (3.3) 

where e is a small constant chosen appropriately. After the evolution of a time interval 
mAt, the orbital point Xi will proceed to ajj+m and neighboring points Xk^ to Xki+m- The 
deviation vectors are thereby mapped to 

yi ^ = 3^fci+m ~ Xi^rnj (3-4) 

which transformation is described by a matrix A\™'\ Then the Lyapunov exponent A is 
calculated by 

A= lim — VlnlU^^ni, (3.5) 

where {n^} is a set of arbitrary unit vectors. If the system is chaotic, A will be positive and 
will not depend on choice of At, e, m and {n*}. 

Fig. 4 shows the Poincare maps for the total angular momentum J = ApM and for several 
values of the spin S = 0.4 ~ lApM by AS" = 0.2pM. As for the spin of a particle with 

mass p, we usually expect that S < 0{p'^). In fact, for a rotating Kerr black hole with mass 
M and angular momentum J, we have the inequality J < M^, where the equality holds 
in the extreme limit. Therefore, for the present case, S/pM = {S/p'^) ■ p/M <, 0{p/M) 
should be much smaller than unity for a test particle, for which we assume that p <C M. 
However, if we would extend our analysis to the case of a relativistic binary system with 
masses mi and m2, we may find that S/pM ~ 0(1) as follows, p and M should be regarded 
as the reduced and total masses, i.e., p = 17111112/ {mi + 777.2), M = f^i + and for the 
spin of the first relativistic star 1, Si/ pM = Si/mim2 = [Si /ml) ■ mi/m2- Then we have 
Si/ pM = Si/m\ < 0(1) for the case of mi = m2. Therefore Si/ pM can be large as unity. 
Because the value of S* = IpM is not a mathematical special bound [Q, we will analyze 
the case of 5* > IpM as well in order to see the spin effect more clearly. 

From Figs. 4a and 4b, which correspond to the type (Bl) potential, we find that no tori 
are broken. Chaos never occurs. On the other hand, in the case of the type (B2) potential, 
which is shown in Figs.4c-4f, we see that some tori are broken, which means that chaos 
occurs for such orbits. This is the first example of chaotic behavior in the motion of a test 
particle in a Schwarzschild spacetime. The sea of chaos spreads in these Poincare sections 
as S increases. We conclude that this chaotic behavior is induced by the spin effect. 

Fig. 5 shows the orbit in the 2-dimensional configuration space corresponding to each 
torus ((1) ~ (3)) in Fig.4d {S = IpM). The orbit (1) is almost perpendicular to the 
equatorial plane. The initial position is near the minimum point of the "effective potential" . 
The orbit (2) is chaotic and the particle approaches the saddle points. The orbit (3) is 
constrained in the very narrow area near the equatorial plane contrary to the orbit (1). The 
initial positions are located near the edge of the "effective potential". In Fig.6, we present 
the Lyapunov exponents A for orbits (1) ~ (3) in Fig. 5. T\/Tp is also given, where Tp is the 
mean orbital period averaged after 10^ rotations around a black hole. T\/Tp gives a naive 
estimation for how many rotations we expect before the chaotic behavior becomes distinct. 
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Wc may justify that the orbit (2) is strongly chaotic because of the large positive vahie of 
its Lyapunov exponent. The reason why the Lyapunov exponents of the orbits (1) and (3) 
converge to positive values may be either that the time series we prepared to calculate the 
Lyapunov exponent arc not sufficient or that the system is not integrable even through we 
cannot see any break down of the tori in the Poincarc map. 

Is there any critical value of the spin for occurrence of chaos? If it is determined by the 
potential type, we obtain Scr — 0.64/xM for J = 4fj,M from Fig. 3. To confirm this, we plot 
the Lyapunov exponent A in terms of the value of the spin S in Fig. 7. After S reaches the 
critical value Scr, A increases remarkably, which supports our conclusion. 

From our numerical investigation, we find several conditions for the occurrence of chaos 
when setting up the initial data. 

• A particle must move in the type (B2) potential. 

• The particle must be bounded in a compact region, i.e., E < Egp, where £^sp = ^(fsp) 
is the potential energy at the saddle point rsp. 

• The particle must have energy E ~ Egp and the position tq, and the angle of the spin 
a must be appropriately chosen in order for the particle to approach the saddle point. 



B. Transition from Regular to Chaos 

Here we shall briefly discuss a mechanism for the occurrence of chaos. 
Fig. 7 shows the Poincare map of five orbits with J = 3.81 fiM, S = IfiM and E — 0.923//. 
Note that E^p = 0.9288/i for this case. We see four different types of orbits. 

• Type (i) : The stable orbits which form a closed curve near the center of the Poincare 
map. 

• Type (ii) : The stable orbits which form a closed curve at the edge of the Poincare 
map. 

• Type (iii) : The stable orbits which form a closed curve in the upper or lower part of 
the Poincare map. 

• Type (iv) : The separatrix of the above three types of orbits. Two hyperbohc fixed 
points in this Poincare section exist. 

It is plausible that the type (iv) orbit plays an important role for the occurrence of chaos. 

The unstable manifold departing from one hyperbolic fixed point will join smoothly to the 
stable manifold of another hyperbolic fixed point. Such an orbit is called a "heteroclinic 
orbit" . 

As E increases from this state, the sea of chaos appears around the hyperbolic fixed 
points and spreads into whole Poincare map. Its behavior is shown in Fig.8. The chaos 
observed here may be caused by a "heteroclinic tangle" , although we need further detailed 
analysis. 

The next question is what causes the existence of type (iv) orbits. Although we do not 
have a definite answer now, we find some correlation between an appearance of such an orbit 
and the type of "effective potential". As S decreases, the type (ii), and then types (iii) and 
(iv) will vanish and only the type (i) orbit remains. This transition is seen in Fig. 4. For the 
case of = 0.8/iM (Fig. 4c), we cannot find the type (ii) orbit, and in the cases of = 0.6 
and OAfiM (Figs. 4a and 4b), all orbits belong to the type (i). The latter case belongs to 
the type (Bl) potential. Hence, it seems that there exists some relation between the type 
of "effective potential" and an appearance of types (n), (iii) and (iv) orbits. We may find 
some criterion for chaos using the "effective potential" . 
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IV. SUMMARY AND DISCUSSION 



In this paper, using the pole-dipole approximation, we study the motion of a spinning 
test particle near a Schwarzschild black hole to clarify its dynamical properties such as a 
chaos. We find the motion of the particle can be chaotic under some appropriate conditions. 
Because the motion of a spinless particle in this spacetime is never chaotic because of its 
integrability, this chaotic behavior is purely induced by the spin-orbit interaction. The 
"effective potential" of the particle is also introduced to classify the dynamical behaviors. 
The "effective potential" s are classified into four different types depending on the total 
angular momentum J and the spin S. When J is large, some orbits are bounded and the 
"effective potential" s are classified into two types: for type (Bl) one saddle point (unstable 
circular orbit) and one minimal point (stable circular orbit) on the equatorial plane exist 
for small spin; and for (B2) two saddle points bifurcate from the equatorial plane and one 
minimal point remains on the equatorial plane for large spin. If J is small, no bound orbits 
exist and the potentials are classified into another two types: for type (Ul) no extremal 
point is found for small spin; and for type (U2) one saddle point appears on the equatorial 
plane, which is unstable in the direction perpendicular to the equatorial plane, for large 
spin. The types (Bl) and (Ul) are the same as those for a spinless particle, but the types 
(B2) and (U2) are new potentials which appear through a spin-orbit coupling. The chaotic 
behavior is found only in the type (B2) potential. We believe that the appearance of saddle 
points is important, because we find chaos only for the orbits which approach the saddle 
points. 

The critical value of the spin beyond which chaos will occur is 5 ~ 0.64/iM for J = 4/iM. 
We also present the Lyapunov exponent, which increases rapidly after the spin 5* gets larger 
than the above critical value. This supports the use of the "effective potential" as a criterion 
for chaos. The typical value of the Lyapunov exponent is about several orbital periods of 
the particle, which may become important in some relativistic astrophysical phenomena. 

In a real astrophysical system such as a binary system, the symmetry of the system is 
lower than the present case. There may be other important effects in addition to the spin- 
orbit interaction, which make the motion more complicated. Then we may expect that chaos 
occurs even in the real system, or that the other effects stabilize the system and chaos will 
never be found. We need further analysis taking into account the other effects such as the 
spin-spin interaction or a force due to multipole moments, even if we adopt a test particle 
analysis |^0|. As for the spin-spin interaction, we can analyze a spinning test particle in 
a rotating Kerr black hole. In our preliminary analysis by the "effective potential" near 
the equatorial plane, the critical value of the spin for chaos gets smaller as the angular 
momentum of the black hole becomes larger. The detailed analysis is under investigation. 

Another important point is whether or not such a chaotic behavior, if it exists, affects 
any realistic astrophysical or physical phenomena. One of the important targets for such an 
investigation is a coalescing binary system, where we need general relativity, in particular 
relativistic dynamics of compact objects as we mentioned in section 1. To examine it, since 
we have studied here only the condition for occurrence of chaos, then we next have to know 
the evolution of the system including emission of gravitational waves. The particle traveling 
around a black hole emits the gravitational waves, extracting the energy and the angular 
momentum from the system. This will tell us whether or not the evolutionary path will get 
into the region of parameter space where chaos will take place. We then have to calculate 

the emission rates of the energy and the angular momentum, E and J for the present system 
[Bl] and follow the evolution. 
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Aid for Scientific Research Fund of the Ministry of Education, Science and Culture (No. 
06302021 and No. 06640412), and by the Waseda University Grant for Special Research 
Projects. 
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Appendix A: The Basic Equations by Use of Spin Vector 



In the text, we have used a spin tensor S'^'^. However, it may be sometimes more 
convenient or more intuitive to describe the basic equations by use of a spin vector S*^, 
which is defined by 



which gives the following constraint: 

P^S, = 0. {A.2) 
The equations of motion ( p.2|) , and ( p.3|) are now 



dr 

Dui" 1 



Dt ^ 
DS" 



where 



The relation between the 4-velocity and the momentum is 

1 

where 

1 1 

* p* — p*a/3 ^ p"/37<5, ( A 

and is the normalization constant determined from v^v'^ = —1. The supplementary 
condition to fix the center of mass is 

VpS^ = 0. {A.9) 

The spin vector is perpendicular to the 4-velocity as we expected. Note that v^S^" ^ 0. 
Equation ( |2.9| ) is just 

52 = S^S^. (A. 10) 

In what follow, we assume a Schwarzschild black hole as the background spacetime. We 
write down explicitly the relation ( |A.1| ) between S'^'^ and S*^ as 
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St 
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r"^ sin e \u' S^''' + u^S'f"' + u'''S' 
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sin 9 
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t ere r nte , 9 ntr 



Conversely, Eqs. (|ATTD ~ (|AT1 ) with Eqs.(|2]|), (|ATID give 

sm t7 



where 



s„ 



Ut 



sin 6' 

Ut 
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Using Eqs.(CT)- (|2^) , Eqs. (|XTT|) ^ (|XT^ ) are now 

1 + fu'^) cos 0" 
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Ug 
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Ug . 
— sm 
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sin 9. 



This expression aheady includes the constants of motion for the anj 

The initial direction of the spin vector at r = rg, 9 = 7i/2, 
equations 



So 
S(h 







Ut 

r 

—1 

Ut 



[J- 



As for the angle of the spin a, we find that 

— 1 Sd 

a = tan 



Using this definition with ([A.IOD , we have 

fiSrl 



Sg 
S^ 



/i^rg + sin^ a 
fiSrI 



: COS a 



/x^rQ + J2 sin^ 



sm a. 



(All) 
(A12) 
{A.13) 
(A. 14) 

(A.15) 
{A.16) 
(A 17) 

(A18) 

(A 19) 
{A.20) 
(A21) 
{A.22) 

;ular momentum. 

^ = is given by the 

(^.23) 

{A.24) 

{A.25) 

{A.26) 

{A.27) 
{A.28) 



a 
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U0,u^ are then some functions of Ut through Eqs. ( |A.24] ) and ([A.25|) . Inserting them into 
the energy conservation equation ( p.24|) , we obtain the quadratic equation for Uf. We also 
have Eq. (|2.8|) for Ur- Solving them, we can finally set up the initial data of and 5*^. 
Setting Pr = Pe = and using Eqs. (|2.8|) and ( |A.10| ), we find the "effective potential" 
. After giving J, S and E and setting initial data, we can solve the dynamical equations 
and ( \A.4\ ) with algebraic equations ( lA.Tp , ( |A.19| ) ~ ( [A.22|) , although here we have 
solved for all variables using the dynamical equations in order to estimate the accuracy by 
the constraint equations. 
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Figure Captions 

FIG. 1: The "effective potential" V(-|-) on the equatorial plane for J = and 5" = 
l/iM. We see that the particle in V(_) moves in the region closer to the event 
horizon (r = 2M) than that in V(+) . 

FIG. 2: Four different types of the "effective potential" V. The saddle point of the 
potential is marked by a cross, (a) Only one saddle point exists on the equa- 
torial plane. This potential is similar to that for a spinless particle. Chaos 
never occur in this case, (b) Two saddle points are found in both sides of the 
equatorial plane. The orbit can be chaotic, (c) When J is very small, the 
centrifugal force is too small to balance with the gravity and then no bound 
region is found. This potential is also similar to that for a spinless particle, 
(d) This also has no bound region. But its shape is different from (c). There 
exists a saddle point on the equatorial plane. This point is locally minimal in 
the r direction but maximal in 9 direction. The particle will eventually fall 
into the black hole after leaving the equatorial plane. 

FIG. 3: The types of the "effective potential" are classified by J and S. The value of 
S at the bottom end of the type (U2) region is slightly smaller than 

FIG. 4: The Poincare Maps for various values of S. All orbits have the total an- 
gular momentum J = 4.0/iM. We set — initially, (a) S — 0.4//M 
and E = 0.97698396/x: The initial position for each torus is Tq = 3.8, 6.0 
and 8.0M. (b) S = 0.6/xM and E = 0.96730999/x and the initial posi- 
tion is ro = 3.6, 5.5, 8.0 and 12.0M. (c) S = 0.8//M, E = 0.95815568// 
and ro = 3.7, 4.2, 6.0 and 8.0M. (d) S = 1.0/xM, E = 0.94738162// and 
(1) ro = 5M, (2) 3.9M(chaotic), (3) 3.72M. The orbits in 2-dimensional 
configuration space are shown Fig.5. (e) S — 1.2/(M, E — 0.93545565/t 
and ro = 3.86, 4.2, 5.2 and 6.0M. (f) S - 1.4//M, E - 0.92292941// and 
ro = 4.5, 5.0, 5.7 and 7.6M. Very strong chaos occurs in this case, although 
it may be unrealistic. 
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FIG. 5: The orbit in the 2-dimensional configuration space corresponding to each torus 
in Fig.4d ((1) ~ (3)). (1) This orbit is almost perpendicular to the equatorial 
plane. The initial position is near the minimum point of the "efi^ective poten- 
tial". (2) The chaotic orbit. The particle approaches the saddle points marked 
by a cross. (3) Contrary to the orbit shown in (1), this orbit is constrained in 
the very narrow area near the equatorial plane. The initial position is located 
near the edge of the "effective potential" . 

FIG. 6: The Lyapunov exponent of each orbit shown in Fig. 5. We find that the orbit 
(2) is strongly chaotic. The ratio of the inverse Lyapunov exponent 1/A to the 
average orbital period t\ = Tp is also shown. 

FIG. 7: The Lyapunov exponent A in terms of the value of spin S. When S becomes 
larger than the critical value S^r ~ 0.64/iM, beyond which we find chaos, A 
increases rapidly. This supports the notion of a critical value of the spin for 
occurrence of chaos. 

FIG. 8: The classification of the orbits for J = 3.81/xM, S = 1/xM and E = 0923/x. We 
find two separatrix (Type (iv)) which divide three types of stable orbits (Type 
(i) ~ (iii))- There exists a heteroclinic orbit starting from one hyperbolic fixed 
point to the other fixed point, which may cause the present chaotic behaviors. 

FIG. 9: When the energy increases from that in Fig.8, the torus gradually spreads out 
around the heteroclinic orbit in the Poincare map. 
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